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TECHNICAL ROTE 2222 


A METHOD FOR THE DETERMINATION OF THE SPANWISE LOAD DISTRIBUTION 
OF A FLEXIBLE SWEPT WING- AT SUBSONIC SPEEDS 
By Richard B. Skoog and Harvey H. Brown 


SUMMARY ‘ 


A method is presented for the determination of the spanwlse load 
distribution of a flexible swept wing at subsonic speeds. The method 
is based on a relaxation approach utilizing aerodynamic loa ding s 
obtained from previously published work (NACA Rep. 921, 19 ^) based on 
Weissinger *s simplified lifting-surface theory together with simple beam 
theory. The solution is expressed in a convenient form such that the 
amount of detailed computing involved when extensive aeroelastic cal- 
culations for many flight conditions are desired is reduced to that for 
a single set of flight conditions. The method is simplified further by 
an abbreviated solution to the relaxation process. Sample computing 
forms and a numerical example are presented to illustrate the method. 


INTRODUCTION 


In the design of unswept wings, the spanwise distribution of .aero- 
dynamic load usually has been considered to be Independent of structur al 
deflections since torsional deflections normally are negligible at 
design operating speeds (usually well below the wing divergence speed) 
and bending deflections are not a factor for zero sweep. On a swept 
wing, however, the span load distribution no longer may be considered to 
be independent of structural deflections since the contribution of 
bending to the streamwise change in section angle of attack can became 
of considerable magnitude as the sweep angle is increased. In addition 
to this aeroelastic effect associated with wing sweep the trend toward 
higher operating dynamic pressures causes factors, which were previously 
negligible, to assume increased significance. 

Several methods have been suggested for calculating the aerodynamic 
loading of the flexible swept wing, both explicitly (references 1 and 2) 
and by iteration (reference 3) • The approach in each of these methods 
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has "been to incorporate aerodynamic and structural theory in an equation 
of equilibrium expressing the balance existing between external and 
internal forces on the wing. To simplify the mathematics involved, 
matrix notations have been employed in two of these methods (references 
1 and 2). Although two of the methods (references 2 and 3 ) indicate how 
aerodynamic theory other than strip theory could be employed, it is 
apparent that such application increases the mathematical difficulties 
a great deal so that the simplest aerodynamic theory has been used in all 
the methods . In addition, although solution of the problem by use of 
matrices has the advantage of simplicity in reducing the necessary com- 
putations to a routine form, the engineer unversed in this mathematical 
tool encounters a loss in physical appreciation. 

The present method arose from the effort to fill the need for a 
mathematically simple approach which can yet include the aerodynamic 
refinements contained in lifting-line or lifting-surface theories. It 
was desired, also, to separate the aeroelastic effects associated with 
the various rigid-win g loadings so that it would not be necessary to 
perform detailed calculations for each new set of flight conditions . 

The method of this report is based on relaxation concepts and employs 
aerodynamic span load distributions from previously published work 
(reference 4) based on Weissinger’s simplified lifting surface theory 
together with structural deflections found from simple beam theory. 

Sample computing forms and calculated effects for an example wing also 
are presented. 

SXMBOIS 


A 

a 


wing aspect 'ratio 

distance from elastic axis to section quarter chord measured 
normal to elastic axis (positive measured forward), inches 
(See fig. 1.) 




ratio of net aerodynamic force along the airplane Z axis 
(positive when directed upward) to the weight of the airplane 

wing span measured normal to plane of symmetry, inches 

section chord parallel to the plane of symmetry, inches 


section chord normal to the elastic axis, inches 

average section chord parallel to the plane of symmetry, inches 


mean aerodynamic chord 


( ^ ) , 

\ I c dy/ 


inches 
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c mo 


c mo 


Cl 

C La 

j 

c l 


section pitching— moment coefficient due to camber as usually 
defined ^section pitching moment ^ 


section pitching-moment coefficient due to camber as used in 
this report (associated with caniber , chord length, and 
dynamic pressure normal to the wing quarter— chord line) 


wing lift coefficient 



rate of change of lift coefficient with angle of attack of 
root section at plane of symmetry 


section lift coefficient 


c l 


a 


section lift coefficient from additional— type loading 



section lift coefficient from basic-type loading associated 
with built-in structural twist 


c i total change in section lift coefficient due to structural 

e deflection 


C l c 


c l 


% 


c l 


e c* 


“o 


C l 


e Az 


changq in section lift coefficient due to structural deflec- 
tion associated with additional, loading 

change in section lift coefficient due to structural deflec- 
tion associated with basic loading 

change in section lift coefficient due to structural deflec- 
tion associated with torsional moment due to camber loading 

change in section lift coefficient due to structural deflec- 
tion associated with inertia loading 


E 


Young 1 s modulus of elasticity, pounds per square inch 




e *b 


total change in section angle of attack due to structural 
deflection, radians 

change in section angle pf attack due to structural deflection 
associated with additional loading, radians 

change in section angle of attack due to structural deflection 
associated with basic loading, radians 



change in section angle of a.ttack due to structural deflection 
associated with torsional moment due to camber loading, radians 
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c hang e in section angle of attack due to structure!, deflection 
associated with inertia loading, radians 

"built-in structural twist of tip section with respect to root 
section, radians 

change in section angle of attack due to structural deflection 
produced by rigid-wing loadings, radians 

c hang e in section angle of attack due to structural deflection 
produced by the aerodynamic loading introduced by deflection, 
radians 

modulus of elasticity in shearing, pounds per square inch. 

distance from elastic axis to section center of gravity 

measured normal to elastic axis (positive measured forward), 
inches 

moment of inertia in bending, inches to the fourth power 
polar moment of inertia, inches to the fourth power 
section lift 'loading, pounds per inch 

section lift loading plus section inertia loading for chord 
sections defined by c* , pounds per inch 

bending moment at any position along the swept span in a 
plane normal to assumed effective root (see fig. l), 
inchr-pound ' 

dimensionless span coordinate [l"/(b/2)J , fraction of semispan 

free— stream, dynamic pressure, pounds per square inch 

wing area (total), square inches 

semispan measured along elastic axis, inches 

torsional, moment at any position al o ng the swept span in a 
plane normal to the elastic axis, inchr-pound 

span coordinate, inches 
(See fig. 1.) 

torsional moment represent ing combined contribution of section 
lift, inertia, anfl moment to the torsional loading about the 
elastic axis at any s panwise station, inch— pounds per inch 
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ang ular change in slope of the elastic axis due to "bending, 
radians 

rotation of wing sections normal, to the elastic azis due to 
torsion, radians 

angle of sweep of the quarter-chord line, degrees 
(See fig. 1.) 

angle of sweep of elastic azis (positive for sweephack), degrees 
(See fig. 1.) 

total airplane weight, pounds 

section structural weight based on chord sections defined by 
c*, pounds per inch 

wing taper ratio 

\root chord/ 


THEORY 

The method discussed in this report is developed to apply to the 
general case of a flexible wing with built-in structural twist, camber, 
and with structural .inertia. The essential feature of the method is 
the application of relaxation 1 procedures to the physical problem of 
dete rmining the aerodynamic span load distribution for the flexible 
wing. In formulat ing the theory in the subsections of the report which 
follow, the unkno wn aerodynamic span load distribution for the flexible 
wing expressed in general terms is applied to the wing together with 
the known load distribution due to inertia. From the bending end tor- 
sional deflections associated with this loading, the rotation (or aero- 
elastic twist) of wing sections parallel to the plane of symmetry then 
is derived. An implicit equation thus is obtained for the twist dis- 
tribution of the flexible wing corresponding to the final, equilibrium 
position ,of the wing- under the combined aerodynamic, elastic, and 
inertia forces. To solve the equation, relaxation methods are applied, 
resulting in a series— type solution for the loading change contributed 
by structural deflection'. As a final step, the lengthy series— type 
solution is converted to gin abbreviated solution which approximates the 
series— type solution very closbly. 

1 In reference 5, relaxation methods are shown to provide a powerful tool 
in solving red un dant problems of structural, frameworks, electrical, 
circuits, and vibratory systems. In the present report the relaxation 
approach employed is based on the principles of that reference rather 
than the exact procedures outlined therein. 


<P 

A c / 4 

Aea 


W 

w 
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In applying relaxation methods to the present problem, the wing is 
assumed to he fixed in position as the rigid-wing loading is applied to 
the wing and then the wing is allowed to deflect under the applied load. 
The wing then is fixed in the deflected position, while the loading 
corresponding to the afore-mentioned deflection as found from rigid-wing 
theory is applied to the wing. The wing then is allowed to deflect 
again. In this way successive deflections can he found which are depend- 
ent on the loading corresponding to the previous deflection. In the 
formulation of the theory, this procedure is used only to develop the 
series— type solution. Actually in practice the abbreviated solution 
previously mentioned is used. In order to apply the method most expedi- 
tiously, available methods for the quick determination of aerodynamic 
span load distributions for wings of arbitrary twist and arbitrary pla n 
form (e.g., reference k) should be used. With any load distribution so 
determined and a knowledge of where to apply the load, it is then a 
simple matter to calculate the amount of structural deflection (bending 
or torsion) due to that load using well-known methods of solution. 

Before discussing the theory in detail, it is desirable.. to establish 
the conventions which are used throughout the report with regard to span 
load distribution. In accordance with the convention of reference 6 the 
spanwise distribution of lift on the rigid wing is considered to be sep- 
arated into an additional lift distribution and a basic lift distribution 
tion, the additional lift being proportional to wing angle of attack and 
the basic lift being dependent only on built— in structural twist. In 
this report, the spanwise lift distribution an the flexible wing is con- 
sidered to be separated into (l) the additional and basic loadings as 
defined above for the rigid wing, and (2) an aeroelastlc loading defined 
as that due to the section angle-of -attack changes produced by structural 
deflection. 

The axes referred to in the report are shown in figure 1. The Y 
axis is the reference axis for all aerodynamic span load distributions 
with the lift assumed to act at the quarter-chord line. The wing is 
assumed to have an effective root perpendicular to the elastic axis and 
located at the intersection of the elastic axis and the plane of sym- 
metry. 

In the following discussion, the basic theory is developed in the 
subsections, loading an Flexible Wing,' Aeroelastlc Twist of Flexible 
Wing, and Evaluation of Aeroelastlc Integrals. The details of this 
development are described in the following steps: 

1. In the subsection. Loading on Flexible Wing, expressions are 
presented for the running load and running torque on the flexible wing 
in terms of the component loading involved. 
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2. In the subsection, Aeroelastic Twist of Flexible Wing, a gen- 
eral expression is developed for the twist due to structural deflection 
based on the unknown aerodynamic loading 2 for the flexible wing and the 
effect of inertia as presented in the first section. It is also shown 
how the general equation may be broken down to lessen the amount of com- 
puting involved when extensive aeroelastic calculations are desired. 

3. In the subsection. Evaluation of Aeroelastic Integrals, it is 
shown how the equation resulting from application of the relaxation 
method can be abbreviated to provide a simple and direct means of evalu- 
ating the implicit twist equations presented in the previous subsection. 

Based on the background developed in these three sections, an equation 
is presented in a fourth subsection. Determination of Span Loading for 
Flexible Wing, showing how the span loading for the flexible wing is 
determined from the various loadings involved. The application of the 
basic theory to determination of lift— curve slope, aerodynamic center, 
and divergence speed for the flexible wing also is discussed in subse- 
quent subsections. 


Loading on Flexible Wing 


The loads which will produce bending of a flexible wing are the 
rigid-wing aerodynamic lift (additional and basic), the inertia load 
normal to the wing, and the increments in aerodynamic lift produced by 
aeroelastic tjrist. The lifting load per unit span based on streamwise 
wing sections can be expressed as 

l - (rigid-wing loading) + (loading produced by wing deflection) 

= ( c i a c <i+ c i- b c <i) + (cz e cq) 

If it is assumed that the effect of taper on sweep is small enough that 
A 0 /4 may be taken equal to A ea and if the effects of inertia or dead 
weight also are included, the lifting load per unit span along the wing 
panel can be expressed as 



2 In setting up the expressions for aeroelastic twist, it is convenient 
to think of the entire flexible wing loading as being unknown. 
Actually all components are known except the component introduced by 
structural deflection. 
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' Ijj- = (c^ cq+ci cq) cos A ea +c^ cq cos A^-wAg (2) 

a d © - 


The aerodynamic and. inertia loadings normal to the plane of the 
win g w ill produce a torsional loading about the elastic axis if their 
lines of application do not coincide with the elastic axis. If a wing 
has camber there will be another torsional loading due to c^. In this ‘ 
report, the contribution of camber has been taken such as to affect only 
torsional deflections of the wing since it was felt that values of 
pitching— moment coefficient for sections normal to the quarter— chord 
line would more likely be" available from two-dimensional tests. In so 
doing, however, it is assumed that simple sweep concepts can be applied 
to finite— span wings without serious error. The torsional moment repre- 
senting the combined contribution of section lift, inertia, and moment 
to the torsional l&ading about the elastic axi^ 3 at any spanwise station 
thus may be expressed as 


tyj- = (rigid-effing torsional loading) + 

(torsional loading produced by wing deflection) 

(c^cq+c^cqja cos A ea -wAgh + c ^ Q c 2 q cob 4 A ea J + 

» 

(c^cqa cos Ae a ) 


Aeroelastic Twist of Flexible Wing 


In thi s section of the report, expressions for the rotation (or 
aeroelastic twist) of streamwise sections of the flexible wing are 
derived in terms of simple integrals of shear and moment based on 


The torsional loading about the elastic axis at tj = rji .(where tjx is 
an arbitrary value of tj) will be a summation of loadings in a pla n e 
normal to the elastic axis. Since this plane will intersect the span- 
wise line of application of the aerodynamic or inertia loadings at 
some point tj = (where i) 2 ^ ijx) an error is introduced, which was 
ignored in the analysis. 

* i 
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elementary beam theory. The general case of the twist due to the load- 
ings on the flexible wing given by equations (2) and (3) is discussed 
first, followed by discussion of a convenient way of breaking down the 
general expression for aeroelastic twist to facilitate computation when 
extensive aeroelastic calculations for a- large number of flight condi- 
tions are desired. • • 


Twist due to loading on flexible wing .— The bending moment at any 
point along the span due to the flexible-wing loading will be 4 


M 


/ *! . P 1 ) 

/ 

.0 J i.o 


drjdrj 


w 


and the torsion will be 


T 


-■ r 


t» in 


( 5 ) 


From elementary beam theory 


11 M A 

* ' s*' 




( 6 ) 


and 




( 7 ) 


For the case of the swept wing, both bending along the elastic axis 
and twisting about the elastic axis cause changes in the streamwise 
angle of attack. The change in section angle of attack 5 due to struc- 
tural deflection for arbitrary bending and twist may be expressed as 


4 Again a slight error is introduced (see footnote 3 ) because the various 
loadings lie alon g different axes. To be correct the loadings should 
be referenced with respect to the elastic axis rather than Y axis. 
The effect of drag on bending moment also was neglected in the anal- 
ysis and normal force was assumed equal to lift force. 

s This involves the assumption of a straight elastic axis for the unde- 
flected wing and ignores the effect of taper in regard to the rotation 
of wing secti ons due to torsion. 



10 


HA.CA. TN 2222 


e = - v sin A ©a + cp cos A ea 


(8) 


In the case of swept-back wings, the two terms, of equation ( 8 ) will "be 
of opposite sign. For swap t-f orward wings, the terms are of like sign. 
Substituting as required in equation ( 8 ), we have 


P 1 ! M pi 1 «n 

e = - s sin A ea J — £ drj + s cos A©a J ^ dT l 


= s 3 sin A 


ea 


/ 


71 -C.o ff.o * di > di > 


El 


dn + 


> ( 9 ) 


r 1 * fjL ** dT1 , 

s 2 cos A ea / “ iT l 

J n GJ 


J 


If substitution is made for l v and t^ from equations (2) and ( 3 ) and 
if the terms are arranged in such a manner as to show clearly the 
various types of loading which contribute to the streamwise twist of a 
flexible wing, the expression for e becomes : 6 


Equation (10) consists of nine terms. The word description given oppo- 
' site each term on page 11 merely gives a physical explanation for the 
existence of the term for the convenience of the reader. 
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. 3 . . rV?Ji.o c la c * ooa Aeefrfa 

e= - s sin A e J — dl l 


_ cj cqa cos A ea dTi 

+ s 2 cos Aea/ — — — § 

G J 


dTi 


0 . rV-Cn-C n c h, cq 003 A ea d n dT l 

- s s sin AeJ — ’ ~ . d n 


J o 


El 


[1] twist due to bend- 
ing caused by 
additional load of 
tbe rigid wing 

[2] twist due to tor- 
sion caused by 
additional load of 
tbe rigid wing 

[3] twist due to bend— 

. Ing caused by 

basic load of 
rigid wing 


+ s 2 cos 




Cj CC[d COS 


J o 


GJ 


dil 


[4J twist due to tor- 
sion caused by 
basic load of 
rigid wing 


+ s 2 


cos 



' ] ^o m o° Sq 0004 AeadT1 
GJ 


dT l 


[5] twist due to tor- 
sion caused by 
cV, of rigid 
wing 


+ s 


3 



wA^dTidtj 

II 


dtj 


[6] twist due to bend— 
ing caused by 
inertia loads 


— s 2 cos 



/i.o vM Z dT l 

GJ 


dr] 


[7] twist due to tor- 
sion caused by 
inertia loads 


— s 3 sin 



[83 twist due to bend- 
ing caused by 
aeroelastic 
loading 


+ s cos 


£ 


/u o c l« cqa cos 


GJ 


dTj 


I 191 
I 

(10) 


twist due to tor- 
sion caused by 
aeroelastic 
loading 


It is evident that tbe first seven terms of tbe prece din g expres- 
sion give tbe twist due to various types of loading of tbe rigid wing . 

Tbe computation of tbis portion of e is fairly straightforward; however. 
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the last two terms (terms [8] and [9]) Include c^ Q (whloh is a f motion 
of e) so that an explicit solution for these terms is not imme diately - 
apparent . 


Twist associated with individual loadings. — As already stated, the 
purpose of the method presented in this report is to determine the aero- 
dynamic spanwise load distribution for a flexible wing. If this load 
distribution is desired for only a single set of flight conditions 
(involving a specific combination of Cl, q, and Jig), the over-all 
effect of flexibility is expressed by equation (10) . However, normally 
it is desired to examine the aeroelastic effects of a wing over a wide 
range of C^, q, and Ag and for various combinations of these factors. 
Therefore, it is usually more convenient to separate c^ Q appearing in 

terms [8] and [9] into the components c, , c, , c, , , and c, 

*e a leb iec^ ^ 

associated with c^ , o^, c^, and A z , respectively. The advantage 

of solving the problem in this way lies in the fact that only one 
detailed computation of the components of c^ 0 need be made. This 
simplification arises from the fact that (l) . c^g is proportional 


to C^, (2) Cj and c^ f do not depend an C^, and (3) 
®b ®mo 



is proportional to A% (or to for a given q and "W/s) . - It is 
possible also to perform the calculations for unit built-in twist and 
unit camber effect. The terms of equation (10) therefore can be 
separated into several components as follows: 


1. Twist for additional loadi ng only.- The aeroelastic twist 
associated with the additional loading, consists of terms [l] 
and [2] of equation (10) plus the portions of terms [8] and [9] 
contributed by c i e& + Since the additional loading is proportional 

to Ql, can be put into the following form: 


ec L = C L[“*aVl A ® a X 


t, f 71 f n C^f—VosAeadTidTi 

'■'x .0 J i.o 


El 


dq + 


qc av 2 s 2 cos A ea / 




GJ 


dq- 


3 r T 

qc av s sin A ea / 

UQ 


^jioQ£) eoBA ^ in + 

EI 


,n f c *e a CN ) 

22 A P 11 ^ 1 »° v^b c av/ 

1 c av s cos A ea / GJ 


cos Aga - — dq 
d q I (11) 


] 
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In this equation the terns within the "brackets are "based on unit C^. 
Once these terms have "been evaluated, Sq. may "be found for any value 
of 0^ simply "by multiplying by 0^. 

2 . Twist for basic loading only. — The aeroelastie twist e e , , 
associated with the basic loading, consists of terms [ 3 ] and [ 4 ] of 
equation ( 10 ) plus the portions of terms [ 8 ] and [9] contributed by 
c l e ^‘ Since the basic loading is proportional to built-in twist, 

e e can be put in the following form: 
b 



+ 


* c av Ss2 cos A ea 


j-tj ^ cos Aq a 7“ dq 

J ^ J i.o\Vay; c av ^ 


GJ 


qc av s 3 sin A ea J 




c 2e. c \ 

) COB ***** 


BI 


: dr] + 


f n f 1 2 p S)-^- cos Aq dTi 
" J x.o VV a ^ Cav 




2 s 2 cos A 




GJ 


>1 


an (12) 


Tn this equation the terms within the brackets are based on unit 
structural twist of the tip section for a given spanwise distribu- 
tion of twist. Once these terms have been evaluated for this dis- 
tribution, e eb can be found for any amount of twist having the 
same distribution simply by* multiplying by e^. 

3, Twist for camber loading nmlv .— The aeroelastie twist ec^, 

associated with the torsional, moment due to cambep loading, consists 
of term [53 of equation {10) plus the portions of terms [8] and [93 
contributed by t . Since the torsion due to camber depends 

upon the amo unt of camber (or the value of cj^),, ®cjjjq can he put 
in the following form, provided the camber is constant across the 
span: 
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r * 

.n Ji 


= Sol {“l c av S s 2 cos Aqq J 


acavs 3 sin J ■ l - tA-La S. 


0.01 1 




COS 4 Ag a dTJ 


dT) — 


COS 'Ag^TjdTl 


dtj + 


qc av 2 s 2 cos 


Q Q 

_/ n 3^- cos AeadTj , 

/ GJ d11 J 

'-'Q 


In this equation the terms within the brackets are based cm a unit 
camber equivalent to c^ = 0.01. Once these terms have been 
evaluated, s c i. may be found for any amount of this type camber 
by multiplyingby cfoo/o.Ol. In the case where a spanwise varia- 
tion in camber is employed, the same procedure adopted for e e-h 
can be used; that is, calculate the terms within the brackets for 
a given camber distribution and then use c^ at some representa- 
tive section as the scale factor. 

4. Twist for inertia loading only .- The aeroelastic twist sa z , 
-i — Ino-rrHn lnnfHricr. consists Of teiHIS [6 J and LT J 


k. Twist for inertia loading only .- The aeroelastic twist sa z . 
associated with the inertia loading, consists of terms [6] and LT J 
of equation (10) plus the portions of terms lo J and 19 J contributed 
bv c l . Since the inertia or dead-weight loading is proportional 
e Ag 

to Ag, 6 a z can be put in the following form: 


e Ag ~ -^Z 


Cav® sin Ag^ 


f- 1 ! -2L- dTjdTi 

■>tl 1.0 i.o Cav 


dn - 


„f <r- <r- d i 

■>11 J i.o av av 


c flT s cos 


qc av s 8 sin Ag a 


p ■! * i »o av av 

W 55-r— a " - 

^aj El 


q c av 2s2 008 ^b. 


r* ff.o C° Ie A Z ) c|? cos Ae ° ai1 . a 1 
■Jo W J 


3 
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In tills equation the terms -within the brackets are based an a 
structural -weight at lg. Once these terms have been evaluated* 

€a 7 may be found far any load factor simply by multiplying by Ag. 
It should be r ememb ered that In combining the effect of Av given 
by thi s equation and the effect of Cj, given by equation (ll), it 
is necessary to adhere to the following relation* in order to retain 
any physical significance: 

. °L4 

A z - WJs 


Evaluation of Aeroelastic Integrals 


Swent-back wings .-- As has already been discussed, terms [8] and 
[9] of equation ( 10) are not immediately solvable; however, as shown in 
the appendix, the solution of these terms can be expressed as a power 
series in q of the following type: 

Ae(n) = fi(Tj)q 2 +f 2(^)9® +f a (n)g 4+ • • • f n(n) 9 n+1 (15) 

Equations (10) through (14) can be expressed, therefore, in series form, 
as : - 


e(ti) = +f 0 h)q+*i(Ti)q 2 +*2(nh s + ^(ti)^ 4 + . . . f n U)q n (16) 


where the coefficients are determined by the particula r requirements of 
the equations (10) through (l4) in mind. In this series (equation (l6) ) 
the values of successive terms are of opposite sign. If equation (l 6) 
is divided by f 0 (t})q, we have 


«<n) . ! + £aM „ + £=<31 + 

f 0(1)9 f o(l) 9 f o(l) 


As has been shown by 0. E. Smith in an unpublished Eorthrup report, 7 
succeeding terms of the series are very nearly related by a constant of 
proportionality so that equation ( 17 ) can be written as: 



1 — kq + k^ 2 . . . 


( 18 ) 


\ similar approach also is contained in outl in e form in reference 7 
under the section titled Effect of Wing Twist. 
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where 8 


fi(n) 

*o(n) 


fg(n) 

fi(n) 


etc. 


As also shown "by Smith, the. series of equation (l8) represents an 
expansion of 1/1+kq so that s(tj) can "be written as: 


e (,) = T*? 

In these equations, f 0 (tj) corresponds to terms [l] through [7] of 
equation (10) (or the first two terms of equation (ll), (12), or (lA-), 
or the first term of equation ( 13 ) ) with q set equal to unity. The 
function f i(tj) corresponds to the twist produced by the aerodynamic 
loading obtained from f 0 (r]) by the method of reference 4-. With f Q (tj) 
.and fi(n) deter min ed, the twist distribution of the flexible wing 
e(rj) can be quickly determined for any value of q by equation (19). 
Then, having e(q), it is a relatively simple matter to get c^ q (tj) by 

the method of reference k, which is baaed on the Weissinger simplified 
lifting-surface theory and. is generalized to permit determination of 
load distribution for a wing of arbitrary plan form and arbitrary con- 
tinuous twist distribution. The reference can be used to provide 
°2 (il) for either a constant lift analysis or a constant angle-of— 
attack analysis. 

It should be noted that the series represented by equation ( 16 ) 
will diverge at same value of dynamic pressure, depending on the struc- 
tural rigidity. An outstanding advantage of equation (19 ) (in addition 
to being brief) Is that no such mathematical difficulty will be encoun- 
tered so that the aeroelastlo effect at any dynamic pressure can be 
calculated. 


8 In practice. It is usually sufficiently accurate to determine k as: 

v - _£±M 

*o(n) 

Since curves of the functions f 0 (n) and fx(Tj) will generally not be 
of exactly the same shape, the value of k will vary somewhat across 
the span. In the case of the example wing, the best approximation was . 
obtained by using the value of k at ti=1.0 since at that station the 
twists given by successive twist distributions (as evaluated using the 
series— type solution) were almost exactly proportional. It should be 
noted that the shape of a given a curve as found from the series-type 
solution will differ slightly, in general, from the shape given by 
f Q (tj) . It lias been found, however, that the differences in curve shape 
encountered do not affect c^ Q (rj) significantly. 
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Swept-f orwar d wings The method can he applied with approximately 
the same accuracy to a swept-? orwar d wing hy rewriting equation (19) so 
that a minus sign appears in the denominator of the function contain- 
ing q. ' 

The equation then becomes 

e (n) = (20) 

the minu s sign arising from the fact that the series for a swept— f orwar d 
wing ( equivalent to equation (18) for a. swept-back wing) is 
l+kq+k 2 q 2 + ... which is merely an expansion of l/l-irq. 


Determination of Aerodynamic Span Loading for Flexible Wing 


The- preceding sections of this report hare laid the background for 
determining the aerodynamic span load distribution for a flexible wing. 
From equation (l) 

i = c la c +. c^c + c 1q O (21) 

From equation (19) it is evident that the distribution of c^ c across 
the span can be written as e 

V 1 ’) (22) 

For q=1.0, this load distribution can be written as 

=T3 f t f o(’l>] (S3) 

© q= i.o 

If solution of f [f 0 (t])] from equation (23) is substituted in equa- 
tion (22) and the resulting expression for c^ Q c is then substituted 
in equation (21), equation (21) can then be written 

2 q(l+k) , . . 

q - °l a ° + "H® + -THET c *e° q=i.o 

If the following relation 

■ c *e = C1 ea + ^ + + ^ 

is substituted in equation (24) and if the various terms are written as 
loading coefficients, equation (24) becomes 
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( 25 ) 


With the various aeroelastic loadings evaluated for q.=1.0, it is a simple 
matter to calculate (l/g.c av )(n) for the flexible wing for any combination, 
of 1* and Az. 


Stability Parameters for Flexible Wing 


T.ift-ciirve sloue .— The lift-curve slope for the flexible wing can 
be found by either of two methods,' depending an whether the viewpoint 
adopted originally was that of constant lift coefficient or constant 
srig ift of attach. Both methods are presented here for convenience: 

1. Constant lift-coefficient analysis .- The aeroelastic span 
load distribution resulting from a constant lift analysis is a 
basic-type loading, which yields zero lift when integrated. In 
solving for the aeroelastic loading by the method of reference 4 , 
the AT»gl «=> of attack of the wing root ay required to obtain zero 
over-all lift also is obtained. This angle of attack represents 
the through which the wing root must be turned in order to 

maintain a given wing lift coefficient for the flexible wing at all 
dynamic pressures. If the lift-curve slope for the rigid wing is 
known, the lift-curve slope for the flexible wing can.be found 
graphically by simply laying off the value of ctj. in the proper 
direction from the rigid-wing lift curve ( increasing wing angle of 
attack for swept— back wings) at the Cj, for which the value of 
was obtained. This proce dur e is indicated in the following sketch 
for swept-back wings: - • 
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The value of is made up of (l) the contrihution associ- 

ated with rigid-wing additional loading and (2) the contribution 
associated with inertia (or dead weigjht) loading, since only the 
contributions associated with these loadings are proportio n a l to 0^. 
These contributions are given by the value b of associated with 

c i e (q) and c^ e (q), respectively. These loadings are found from 
a Ag 

equations (ll) and (l4) previously presented. Lift-curve slope for 
rigid wings is given directly in reference 4 for a wide range of 
plan forms . 


2 . Constant angle— of-att ac k analysis .— The aeroelastic span 
load distribution resulting from a constant angle— of— attack analysis 
is similar to an additional-type loading and yields a lift when 
integrated. For a swept-back wing,, the c h a ng e in lift resulting from 
integration of the aeroelastic loading associated with the twist of 
equation ( 19 ) can. be written as 


dC L = 



(2 6 ) 


where 

dCiQ lift coefficient resulting from integration of the aeroelastic 
loading corresponding to e(q) 

£Cjj X lift coefficient resulting from integration of the aeroelastic 
loading corresponding to f 0 (q) 
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The lift coefficient for the flexible wing ( swept-back) therefore 
is 

^ - xSe ' < 2 t> 


where 


Cl 0 lift coefficient of the rigid wing 

If equation (27) is divided by the arbitrary angle of attack, the 
lift-curve slope for the flexible wing can be written as 


(Cl teV (Cl to R C 1 1«Sl C^) 


(28) 


where the subscripts F and R refer to the flexible and rigid 
wing, respectively, in the case of a swept-forward wing, equation 
(28) becomes 


(GL a } F ^^R ( 1+ 1-ka c/) 


(29) 


Aerodynamic center .— The accepted definition of aerodynamic center 
of a rigid wing is the centroid of all the additional loads. It is evi- 
dent that an a flexible wing the lift increment due to angle of. attack 
includes not only an additional-type load (in the rigid-wing sense) but 
also a varying amount of aeroelastic lift. The effective aerodynamic 
center of a flexible wing will thus differ from that of the rigid wing 
but will still be the centroid of 'all the additional loads . In stability 
analyses, it usually is customary to neglect the vertical location of the 
centroid since the effect of drag on stability normally is negligible. 


The varying amount of aeroelastic lift is made up of (l) the com- 
ponent associated with, the additional loading and (2) the component 
associated with the inertia (or dead weight) loading, since only the 
components associated with these loadings are proportional to C^. 

These components, c 7 (n) a n d oi (tj), respectively, are defined by 

®a ®A2 

the twist distributions given by equations (ll) and (lk) previously pre- 
sented. With Ci (n) and ci (ij) determined, the aerody nami c can— 
^a , . e A z 

ter of the flexible wing can be found by any method for determining 
centroids, remembering that the chardwise load is assumed to lie along 
the quarter— chord line of the wing. The same method is applicable, in 
general, whether the analysis adopted is for cons tant lif t or constant 
angle of attack. 
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Divergence Speed for Flexible Bing 


The ease with, which divergence speed can be obtained from equations 
(19) and (20) is worthy of brief mention also. If equation (19) is dif- 
ferentiated with respect to q, the derivative becomes 

k = (55JF ^ l30> 

The condition for divergence is that the derivative (d/dq) e( tj) must 
equal infinity. To satisfy this condition. 


1 + kq D = 0 (31) 

where qjj is the divergence dynamic pressure, so that 

< 3 d = “ ie (32) 

In the case of a swept-forward wing, differentiating equation (20) and 
proceeding as before yields 

«d - ; (33) 

Equations (32) and (33) show the familiar fact that a swepb-forward wing 
will diverge at son© positive value of q while a swept— back wing will 
not diverge at any positive value of q. 


AEELICATIOB 


Computing Forms 


Computing forms to aid in solving for the change in span load dis- 
tribution due to wing flexibility are presented in table I. The basic 
data needed to perform the calculations are shown in table I (a). Tables 
1(b), 1(c), 1(d), and 1(e) are far the purpose of computing the functions 
(f 0 (tl) and f x ( tj ) ) for each of the rigid-^ring loadings, aid require essen- 
tially the same operations. Any differences noted are merely for ease in 
Tm-nfl 1 -tng the c ompu tations for a given load i ng. Table l(f) is essentially 
the same as table l(c) which is for basic-type loadings. The only cLis r 
tinction between the two forms is that table 1(f) is far the aeroelastic 
loading introduced by wing deflection and table l(c) is far the basic 
loading due to built-in twist. With one exception these farms consist of 
separate calculations of bending deflection (columns 4 through 7) s^^d of 
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torsional deflection (columns 9 through. 11) . ' 'ghe single exception iB 
table l(e) for camber loading which, does not require bending com- 

putations for obvious reasons. The last column (column 12) in each 
case is the sum of the bending and torsional deflections in terms of 
streanwise angl e — of -attack change. This column yields f q (t|) in the 
case of tables l(b), l(c), l(d), and l(e) and yields fi(Tj) in the 
case of table l(f). The column headings are either self-explanatory or 
are explained in the computing instructions following the tables. In 
the computing instructions operation. A integrates the r unning load 
normal to the wing to obtain the shear at designated spanwise stations. 
Operation B Integrates the shear curve so obtained to dete rmine the 
bending moment at the same spanwise stations. Operation D integrates 
the running torsional load to obt ain the torsional moment at the same 
spanwise stations. The summations are performed as indicated from 
H = 1.0 to t) = 0 . Operation C integrates the M/El curve to obtain 
the slope of the elastic axis at the chosen spanwise stations. Opera- 
tion E integrates the T/GJ curve to obtain the twist about the elastic 
axis at the same spanwise stations. The summations in these two opera- 
tions are performed as indicated from tj = 0 to r\ = 1.0. The 
integrating operations have been set up in accordance with the trapezoi— , 
dal rule for approximate integration. The spanwise stations used in 
the computations, therefore, should be of sufficient number and of 
adequate spacing so that the integrations will not be subject to signif- 
icant error. 


numerical Example 


The solution procedure Indicated in the preceding section of this 
’report has been applied to a relatively flexible example wing of moder- 
ate sweep and high aspect ratio. Compressibility considerations were 
neglected in the present example, since the modifying influence of com- 
pressibility is small compared to the isolated, primary influence of 
dynamic pressure. The geometry of the example wing is shown in figure 2 
together with curves of the structural— stiffness data used in the cal- 
culations. The wing has an aspect ratio of 9*43, a taper ratio of 0.42, 
and a sweep angle of the quarter-chord line of 35° • As indicated in 
the figure, the elastic axis is located at 3&-percent chord. 

The rigid-wing loading curves used in the calculations are shown in 
figure 3. The aerodynamic loadings (additional, and basic) were obtained 
from reference 4, neglecting compressibility. 9 The additio n al— type 

9 "" — — — — ■ — - 1 

Eor the reader interested in including compressibility, it should be 
noted that the method of reference 4 accounts for compressibility 
effects on the basis of linearized compressible flow theory so that 
the modifying influence of compressibility can easily be included 
using that reference by merely following the procedure outlined 
therein. 
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loading corresponds to an over-ell wing lift coefficient of 1.0 and the 
basic loading corresponds to a linear twist distribution hawing 1.0° 
washout at the wing tip. The dead-weight shear distribution was obtained- 
from weight estimates for the wing. The abrupt shear changes shown are 
due to the concentrated masses of the engines. The camber loading was 
obtained assuming constant camber across the span equivalent to 
CjjJq — 0.01. 

Loading. — Computations for the example wing assuming only the addi- 
tional loading to exist and based on a oonstant lift analysis are pre- 
sented in table II to illustrate use of the computing forms. Table 
H(a) presents the geometric and structural parameters for the example 
wing for specified stations along the semispan. Table H(b) presents 
the calculations for the angle-of -attack redistribution from which the 
function f 0 (q) is found. The function fi(ti) is found from the 
angle-of-attack redistribution calculated in table Il(c) . In table 
11(b), the loading used is the additional loading. In table H(c), the 
loading used is the basic loading found from the twist distribution cal- 
culated in table Il(b). The functions f 0 (ij) and f ( t|) and the ratio 
are Plotted in figure h against spanwise station f[. The 
ratio fi(ti)/fo(Ti) Is plotted for the purpose of illustrating the 
differences in curve shape between the functions f 0 (Tj) and fi(ri). 

As can be seen from the figure, the distributions of twist (as defined 
by the shape of the curves) for f 0 (^) and fi(q) are somewhat dif- 
ferent due to the fact that f 0 (n) was determined from an additional 
loading and fi(ri) was determined from a "basic loading. In spite of 
the difference shown, however, the assumption of proportionality between 
successive terms of the series is sufficiently close to provide a good 
approximation since basic— type loading is affected very little hy dif- 
ferences in curve shape such as shown in the figure. As can he seen 
from the figure, a large variation. in fi(ii)/f 0 (ii) " across the semispan 
can exist so that it is important to choose the value of this ratio at 
the proper value of tj. As stated earlier, the best approximation to 
the series— type solution, was obtained for the example wing by choosing 
the value of fi(Ti)/f 0 (il) at the tip station. These remarks also apply 
if a constant angle-of-attack analysis is adopted. 

The aeroelastic loadings associated with each of the rigid-wing 
loadings are presented in figure 5 for several values of dynamic pres- 
sure as obtained from the constant— lift analysis. With the solution in 
the form shown in the figure, the computation of the spanwise load 
distribution for the flexible wing can be found relatively simple for 
any set of flight conditions and any set of camber and twist distribu- 
tion similar to that assumed merely hy combining the loading for the 
rigid wing with the proper combination of basic loadings due to deflec- 
tion as indicated by -equation (25). As has already been pointed out, 
such a computing shortcut is made possible by the linear aerodynamic 
and structural theory of the analysis which renders the deflection 
loadings (at a given dynamic pressure) proportional to either 0^, A & 
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e^, or c^, depending on the rigid-wing loading involved. The aero- 
elastic loadings for a constant angle— of— attack analysis are not pre- 
sented since the proper combination of aeroelastic loadings is believed 
to be indicated sufficiently by figure 5. 

Stability parameters .— Since no experimental data were available 
with which theoretical span load distributions of a flexible wing could 
be compared, the validity of the method was checked to some extent by 
comparing predictions of lift— curve— slope change and aerodynamic— center 
shift with those predicted by the method of reference 1. 

The variation in lift-curve slope with dynamic pressure for the 
example wing, as calculated by the procedures outlined earlier in the 
report, is presented in figure 6, neglecting the modifying influence of 
compressibility. In the calculations, a rigid-wing lift-curve slope of 
0.071 obtained from reference 4 was used. In the figure, a curve 
showing the slope variation according to the method of reference 1 also 
is presented for comparison. As can be seen, agreement between the two 
methods is very good. 

The spanwise shift in the centroid of the additional loads with . 
dynamic pressure for one panel of the example wing is presented in 
figure 7 together with the corresponding shift in aerodynamic center 
parallel to the plane of symmetry. At a dynamic pressure of 500 pounds 
per square foot, the spanwise shift is shown to be inboard by about 
6 percent of the semispan. The corresponding chordwise shift is shown 
to be forward by about 20 percent of the mean aerodynamic chord. Good 
agreement with the prediction of reference 1 is shown again. 


Ames' Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Calif., Judy 31, 1950. 
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APPENDIX 

DEKEVATIOn OF WING TWIST POWER SEKEES 

In the text of this report, the streamwise twist of a flexible, 
swept wing is given by equation (10) as a summation of terms due to the 
various rigid-wing loadings .and the aeroelastic loading introduced by 
flexibility. Equation ( 10 ) can be summarized as follows: 

6 (ti) = e 0 (q)-tAe (n) (Al) 

where 

e 0 (q) the twist of the flexible wing due to the total rigi during 

loading (as given by terms [l] through [7] of equation (10)) 

and 

Ae(q) the twist of the flexible wing due to the aeroelastic loading (as 
given by terms [ 8 ] and [9! of equation ( 10 )) 

As stated in the text, the value of Ae(q) in the above equation 
cannot be evaluated directly In any explicit manner so that a method of 
successive approximations must be employed. 

One method for evaluating e(tj) is to adopt the obvious iterative 
approach and compute e(q) by successive approximations to Ae(q) until 
sufficient accuracy is obtained. A more convenient method, however, is 
to apply a relaxation procedure to the determination of Ae(q) . In this 
method, the wing is assumed to be fixed in position as the rigid-^wing 
loading is applied, and then the wing is allowed to deflect under the 
applied load, resulting in the streamwise twist distribution e^q) of 
equation (Al) . The wing then is fixed in position again and the rigid- 
wing loading is removed. The aerodynamic loading corresponding to 
6 Q (q) then is applied to the wing and again the wing is allowed to 
defleot (in accordance with the applied load only), resulti n g in a new 
twist distribution. In this way successive twist distributions can be 
found which are dependent upon the aerodynamic loading corresponding to 
the previous twist distribution. The final twist distribution for the 
flexible wing can therefore be expressed as 

6 (tj) = e 0 (q)-fAei(q)-4Ae 2 (q) . . . + Ae^r^) (a 2 ) 


Comparison of equation (A 2 ) with equation (Al) shows that 
Ae(q) = A€ 1 ( n ).-tAe 2 (q) . . . + A^(q) 


(A3) 
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New since® 


where 


Aeih) = ql 


c Ze ^(n) = f ^oCnij 


GnU) = 11 ’ 


|^c z (ti)J 


€ o(n) = <i Mn) 

it is apparent that equation (A4) can he written as 

Ae x (Ti) * q^fi(q) 

Similarly, 

. A^(tj) = q I-j" c 2e2 (Tl)J 

where 

c l ez W = f [^id)] 

so that, with equation (A5), it is apparent that 


In like manner 


A€ 2 (ti) = q f 2 (n) 


Ae n (ti) - q n+1 f n (t]) 


9 In these expressions, the notations f £ J • and J merely 

indi cate the general dependence of loading and twist on the associated 


twist and loading, respectively. The notation f- 


a specific function. 


n 


indicates 
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It should he noted that for a swept-back wing each new twist dis- 
tribution calculation by this procedure will be of opposite sign to the 
twist immediately preceding, since successive aerodynamic loadings will 
exert bending moments of apposite sense due to the bending-twist rela- 
tionship for a swept— back wing. Summing the successive contributions 
to Ae(ti) of equation (Al)- as given by equations (A5)> (A7)> and (A8), 
we have 

Ae(rj) = q2fi(rj) '+ q 3 f 2 (Tl) + . . . q. n+:i Tja(n) (A9) 

Substituting in equation (Al) for e 0 (Tj) by the relation given under 
equation (A?) and for A€(tj) as given by equation (A9), it can be .seen 
that the twist distribution for the flexible wing can be expressed as a 
power series in q as follows: 


e (n) = q.f 0 (n) + AlU) + i 3 f 2 (n) + . . . ^^(n) (aio; 
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TAHU5 I.— COHTHTOED 
(a.) Dead weight loading 
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(e) Camber loading 
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03ABLE H.- SAMPLE COMFDTATICBT FCR EXAMPLE WING BASED OH ADDITICHAI/-TIEI5 LOADIHG 

4 

(a) Geonetrlo and structural parameters. 
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(fc). Calculation of the function f 0 (t|) 
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TABUS II.- COMCLUOED 
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(c) Calculation of tiue function fi(q). 
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Figure /.- Sketch of swept -back wing sho 
in analysis. 
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Figure 2 Sketch of example wing with pertinent 
dimensions, plan-form parameters, and elastic 
properties. 




Spanwise station , y 
(d) Pitching moment di/e to camber. 

Figure 3.- Rigid -wing loadings fdr the^ example Wing. - 
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Figure 4.- The twist distributions fo(ip) and ft (ip) and 

ft (if) 

the ratio — — for the example wing as caicu- 
fc(v) 

iated for the additional loading . 
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. (a) Additional loading. 

Figure 5.- Aeroeiastic loadings and deflections associ- 
ated with the various rigid -wing loadings for the 
example wing as calculated for several values of 
dynamic pressure using a constant -lift analysis. 
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(c) Dead -weight shear. 


Figure 5- Continued. 
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Figure 7.- Variation in spanwise and chord wise positions of aerodynamic center 
with dynamic pressure for example wing. 
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